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We study the magnetization plateau state of the three- leg spin-| tube in the strong rung coupling 
region, where S'3-symmetry breakings and the low-energy chirality degree of freedom play crucial 
roles. On the basis of the effective chirality model and density matrix renormalization group, we 
clarify that, as the leg coupling increases, the chirality liquid with gapless non-magnetic excitations, 
the spin imbalance phase and the vector-spin-chirality ordered phase emerge without closing the 
plateau spin gap. The relevance of these results to experiments is also discussed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Kz, 75.40.Cx 
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I. INTRODUCTION 

Geometrical frustration on magnetism has long been 
one of the attractive subjects in condensed-matter 
and statistical physics, since the frustration provides 
rich physical phenomena and various ordered/disordered 
state&i'. It is well-established that the spin chirality of- 
ten plays a fundamental role as we probe the frustration 
effects, especially, in the triangular lattice systems^ii^. 
Recently multiple-spin orders without any magnetic mo- 
ment, including vector spin chiral order, have been ac- 
tively studied as a new topic in frustrated magnetism 
(e.g., one- and two-dimensional J1-J2 spin models^"^). 
The vector spin chirality also attracts extensive attention 
in the context of multiferroics^, where the chirality order 
induces electric polarization. In the most of frustrated 
systems like the J1-J2 models, however, the chirality ex- 
citation is usually embedded in conventional magnetic 
excitations, which make direct observation of the chiral- 
ity difficult. In order to gain deeper understanding of the 
frustration physics, thus, it may be a key issue to extract 
the chirality excitation energetically separated from the 
magnetic fluctuations in a realistic situation. 

Among a mount of frustrating systems, the three- 
leg spin tube, consisting of coupled three spin-i an- 
tiferromagnetic chains [see Fig. [T] (a)], is one of the 
models deeply related to spin chirality; We can de- 
fine clockwise/anticlockwise rotation along the rung in 
the spin tube. In fact, the topological structure of 
the spin tube is known to induce several interesting 
phenomena^Ti^^. Recently, spin-tube materials such as 
[(CuCl2tachH)3Cl]Cl22ii22 CsCrF^^^iSi have been re- 
ally synthesized and characteristic properties to the spin 
tube have been revealed by several experimental ap- 
proaches. In particular, it is pointed out that the broad 
peak of specific heat is associated with a gapful chirality 
excitation in the twisted tube [(CuCl2tachII)3Cl]Cl2^. 
However, it should be also noted that the contribu- 
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FIG. 1. (color online) (a) Structure of the three-leg spin tube 
and (b) ground-state phase diagram of the i plateau state of 
the spin tube. The plateau is predicted to vanish at a strong 
leg-coupling point J/ J' ~ 0. 1^^'^'^ (see the text). 



tion from gapless magnetic excitation overlaps this broad 
peak related to chirality. 

In this paper, we demonstrate that the quantum phase 
transitions associated with the chirality actually occur in 
the magnetization plateau of the straight quantum spin 
tube, where energy scale of the chirality is certainly sepa- 
rated from gapful magnetic excitations. The Hamiltonian 
of the spin tube is given by 
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the spin-i matrix, J(J') > is the 
intra(inter)-triangle coupling, and i (j) represents the la- 
bel of the rung (leg) direction (i: mod 3). This model 
(H)) looks very simple, but the frustration due to the tube 
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structure is expected to induce various characteristic 
properties. In fact, it was shown that the model ([1]) has 
a uniform vector spin chirahty order in the weak rung- 
couphng region (J ^ J') in a magnetic field H^^-^^. A 
rather interesting parameter region is the strong-coupling 
limit (J ^ J'), where the system is basically described 
by the weakly coupled triangles. In the strong rung limit, 
the composite spin 
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on each unit triangle is classified into T ~ | © ^ © ^ 
sectors and then the = i states of T = i sectors 
lead to a robust magnetization plateau at i of the full 
moment. A key point is that the two-fold degeneracy 
of T = i sectors in this plateau state brings an active 
low-energy variable, which is just the chirality degree of 
freedom. Utilizing the low-energy effective model and 
density matrix renormalization group (DMRG), we will 
show that the energetic separation of the spin and chiral- 
ity excitations leads to nontrivial quantum phase transi- 
tions without destroying the magnetization plateau. The 
main results are summarized in Fig.[l](b); we find chiral- 
ity liquid, spin imbalance, and the ferro-chirality ordered 
phases. We also explain that these orders are accompa- 
nied by the S'3-symmetry breaking in the unit triangle. 

The remaining part of this paper is organized as fol- 
lows. In SecHH we study the ^ plateau state based on the 
effective spin chirality model. We also discuss the role of 
the S'3-symmetry in the quantum spin tube. Section IIIII 
is devoted to the numerical results derived from DMRG 
method. Combining the DMRG results with the analyt- 
ical predictions in Sec. |lTl we reveal three new phases in 
the plateau region; chirality liquid, spin-imbalance, and 
the ferro-chirality ordered phases. Finally we summarize 
our result and the relation between it and previous stud- 
ies in Sec. IIVI Furthermore, we discuss the relevance of 
our result to experiments. 



II. EFFECTIVE CHIRALITY MODEL AND 5^3 
SYMMETRY 



Let us start with the low-energy effective theory for the 
plateau state in the strong rung-coupling region. We can 
represent the two-fold degenerating bases for the — 
1/2 states of T = i on each triangle as 
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g27rj/3 ^ ^^^^ denotes the left- (right- 
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handed mode in the rung direction^. These two states 
indeed stand for the chirality degree of freedom. By pro- 
jecting out the high energy states with — and 
T = I in every unit triangle, the effective Hamiltonian 



of the plateau state is obtained as 
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(2) where Tj is the pseudo-spin-^ matrix defined by tJ = 
{\L)jj {L\ ~ \R)jj {R\)/2. The coupling constants are eval- 
uated as K^y = 2J'/3 ~ 5J'V(27J), ^ -J'^/J, 
K'^y = 8J'V(27J) and ATg = -16J'V(27J) within the 
second-order perturbation in J' . Here, it is worthy to 



note that the relation between Tj and Si,j is given by 
— \/iPjXjPj and — —PjfijPj, where 
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are respectively the z component of the vector spin chi- 
rality and an imbalanced magnetization on each triangle, 
and Pj = \L)jj{L\ + \R)jj{R\ is the projection operator 
to the T/ = i states of T = i. 

In order to resolve possible quantum phase transitions, 
it is very instructive to discuss the discrete symmetry of 
the spin tube. The spin tube has iSa-group (= Csv point 
group) symmetry in the rung direction in addition to the 
translational symmetry along the leg direction. The op- 
erations in the S3 group are composed of the cyclic per- 



mutation Si 



with mod 3 and the reflection 



Si,j -f-?- Si'j at a bond in every unit triangle (i ^ i'). 
Possible iSs-symmetry breakings are classified by its sub- 
groups: (a) the bond-parity breaking with conserving the 
cyclic symmetry, (b) the cyclic Z3 symmetry breaking 
with conserving a part of bond-party symmetry, or (c) 
the full breaking of the S3 symmetry. The vector spin 
chirality Xj is a typical order parameter in the case (a), 
which changes its sign by the reflection, but is invariant 
under the cyclic permutation. This cyclic symmetry is re- 
lated to the spin current circulating in the rung direction. 
On the other hand, /ij can be an order parameter of the 
case (b), since its form changes via the cyclic permuta- 
tion, but is invariant under the reflection S2,j ^ ^Sj- If 
II j becomes flnite, it suggests that the isosceles-triangle- 
type imbalance occurs for {Sf j) in the plateau state. 

We discuss the relation between the 5*3 symmetry and 
the effective model (|H) . Write the cyclic permutation op- 
eration of the S3 symmetry group as Tc, and the bond 
reflection as 7^(= 7^^^). In the level of the effective chi- 
rality T, the S3 symmetric operation is given by 



TrT+Tr = 



n- +n — 1 2 + 



TrT^ Tr = T+ 



(6) 



for any j. Under these operations of the S3 symmetry, 
the effective Hamiltonian Q is confirmed to be invariant. 
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Here we should remark that in the model (jj]) , the second- 
roder perturbation process generates the U(l)-symmetry 
breaking term, although the U(l)-symmetric XY 
model, which is obtained within the first-order pertur- 
bation, has been often used for the spin tubesiSiiiii^. 
This is consistent with the fact that rj ~ Xj is not ex- 
actly conserved in the original spin tube. Thus we need 
a careful consideration about the role of symmetry and 
interactions in the effective model (|4|). 

According to the bosonization approach^^, the low- 
energy physics of the model (HJ is described by a massless 
free boson theory with several interactions. The effec- 
tive Hamiltonian for the free boson, i.e., the Tomonaga- 
Luttinger (TL) liquid is represented as 



can have a finite expectation value through the relation 



"Htl = / dx- 



k{djf + k-\d.. 
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where (</>, 9) is the canonical pair of scalar fields [x — ja 
and a is lattice spacing), K is the TL-liquid parame- 
ter, and V is the low-energy excitation velocity of the 
model (U]). The effective spin tj and the bosonic fields 
(6, 9) is related as 



a d(f){x) 
x/tt dx 
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with non- universal constants ai, bo and bi. The 6*3- 
symmetry operations on the effective fields are summa- 
rized as 

Tre{x)Tr = ~9{x), Tr(l){x)Tr = ~4>{x) + 0F/2, 

Tc9{x)T~^ = e{x) + 27^/3. (9) 

In addition, the operation of one-site translation along 
the leg 71 transforms the boson fields as 



Tie{x)Tf^ = 6l(x + a) + v^, 

ri(f>{x)rr'^ = <j){x + a) + 



(10) 



These symmetries impose significant restriction to the 
possible interaction terms in the effective field theory. 
Among various vertex operators permitted by the ^3 and 
translational symmetries, the most relevant terms are 
given by for cos(2-\/27r(/)) and cos{6^/n6), for which the 
scaling dimensions are respectively AK and 9/K. Since 
the value of K approaches unity in the J' / J — limit 
(the XY model) , we can see that the interaction terms in 
Eq. (U]) are all irrelevant for sufficiently small J', suggest- 
ing that the critical chirality liquid is realized in a certain 
region of small J'. On the other hand, the system may 
have two kind of instabilities as J' increases. The first 
case is the ferro-chirality order of rj ^ Xj- Since the 
negative in Eq. ^ raises the value of K to +00, 
the ferromagnetic instability may occur, at which the ve- 
locity V also vanishes. The other case is the staggered 
order of the imbalanced magnetization /z^; If 9/K < 2, 
6'-field is locked and then the staggered component of 



{—y cos{y/7T9). Here, we note that, in 



the following numerical computations, the ferro-chirality 
oder actually appears, but a uniform order of /ij is real- 
ized rather than the staggered type. 



III. NUMERICAL RESULTS 

Now we apply DMRG to the spin tube model ([1]) to 
quantitatively examine the transitions and orderings with 
the help of results in SecHH We fix J = 1 in the following 
numerical calculations. 



A. chirality liquid phase 

First, we focus on a sufficiently strong-rung coupling 
region. In Fig. [21 we present the longitudinal spin cor- 
relation function {S^jS^j,) for L = 96 systems with 
J' = 0.01, • • • , 0.45. The rapid decay near the right edge 
in Fig. [2] comes from the open boundary effect. Thus 
it can be confirmed that the correlation function follows 
a power-law decay for \j — j'| < 50: {S^^Sfy) — m? ^ 

{—y^-' \j — , where m = ^ is the uniform magneti- 
zation per spin and 77 is the critical exponent. This decay 
fashion is in agreement with the prediction from the effec- 
tive TL-liquid theory ([7]). We can also see that 77 becomes 
close to 0.5 in the J' — limit, where the Hamiltonian 
reduces to the XY model. As J' increases, 77 approaches 
zero toward the ferro-chirality transition. Utilizing the 
effective field theory ([7]) based on the XXZ chain we 
can evaluate the critical exponent rj in the strong rung- 
coupling region J ^ J' . The value upto the second or- 
der of J' is given by r/ ~ 0.5 - 0.885 J' -I- 0.640 J'^ + ■■■, 
where we have assumed the nonuniversal parameter &o — 
0.5424 • • •— . We have confirmed that this value of r/ is 
semi-quantitatively consistent with the numerically esti- 
mated value from the correlation function of Fig. [5] in 
J ^ J'. From these results, we conclude that the gap- 
less non-magnetic chirality excitation is described by the 
effective model (|4]). Here, note that the width of the 
plateau is sufficiently large for J' < 0.5 and the trans- 
verse correlator (SfjSfj,) exponentially decays, indicat- 
ing that the magnetic excitation has a large gap corre- 
sponding to the plateau width. 



B. ordered phases 

As J' further increases, the negative derives the 
system toward a ferro-chirality ordered state with (xj) ^ 
0. Figure [3] illustrates the results of the order parameters 
X — iXj) and = {^ij). Here, x is observed at the center 
triangle of the tube of size L = 120^ and ^ is the bulk ex- 
pectation value based on the infinite system DMRG. We 
have checked that the boundary effect is negligible within 
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FIG. 3. (color online) Expectation value of the order param- 
eters X (solid circle) and /i (open circle). Inset shows these 

FIG. 2. (color online) Longitudinal spin correlation function ^^^^^ parameters around the transition points. The vertical 

i^lj^j') of the spin tube for J' = 0.01, 0.1, 0.2, 0.3, 0.4, broken lines indicate the transition points J^i and 

and 0.45 from bottom to top, where m = |. Two sohd lines 

indicate guides for 77 = 0.5 (XY chain) and 0.1. 



computations for L = 96, 120 and 144. From the main 
panel, we can see two quantum phase transitions near 
J' = 0.5. Note that the plateau width around J' = 0.5 
is about 0.5 J, which is sufficiently larger than the energy 
scale of the non-magnetic chirality excitation. Figure [3] 
clearly shows the emergence of the ferro-chirality order in 
J' > •^c2 = 0.496, which is consistent with the effective 
model (jlf. We have confirmed that this ferro-chirality 
order extends to J' > 1 and thus it would be adiabati- 
cally connected to the vector chirality order in the region 
of the weakly-coupled three chainsi^. Here, note that 
both {SfjSfj,) and {S^^Sf^,) show exponential decays 
in J' > J^2 thus the magnetic and chirality excita- 
tions have finite gaps in this chirality ordered phase. 

From the inset of Fig. |3l we also find that the spin im- 
balance phase emerges in a narrow region J'^^ < J' < 
with J'^^ ~ 0.478. In this region, the symmetry of the unit 
triangle reduces to the isosceles type, where the expecta- 
tion value of one spin of each rung triangle is larger than 
those of the remaining two spins: {S^ j) > {Sf_^_i j) = 
{^i+2,j)- Fig. m we present the (Sfj) distribution for 
J' = 0.485, which exhibits a typical spin profile of the 
spin-imbalance state. The open-boundary effect rapidly 
decays and a uniform spin imbalance along the chain di- 
rection is realized around the center of the tube. Figure[5] 
shows a semi-log plot of (SfjS^j,) — mf, where rui is the 
bulk expectation value of S'f ^ calculated at the center of 
the tube. The exponential decay of the correlation func- 
tions in Fig. [5] indicates that the system is gapful. We 
note that the imbalanced nature is present not only in 
the magnetization profile, but also in the spin correla- 
tion functions. As we see from the inset of Fig. [SJ the 
correlation length for the less polarized spins becomes 
divergent as J' J'^i + 0, while that for the most po- 
larized spin remains finite value. This suggests that the 
instability of the spin imbalance toward the chirality liq- 



uid state (J' < J'^i) may be governed by the fluctuation 
of the less polarized spins of the triangle, although the 
critical behavior of fi cannot be determined within the 
accuracy of the present DMRG results. As J' increases, 
the correlation lengths of the most polarized spin and the 
remaining two become comparable with each other and 
finally arrives at the ferro-chirality transition point Jj^,. 
Here, we note that, for 0.485 ^ J' < Jc2i the spin correla- 
tion functions becomes highly oscillating and thus precise 
estimation of the correlation length is difficult. We stress 
that this imbalanced order cannot be described by the ef- 
fective model (|4]). This suggests that the hybridization of 
= 3/2 sector plays an essential role in the imbalanced 
phase (see the following paragraphs). On the other hand, 
the jump of the order parameters at clearly shows 
that the transition at J' = of ^''^t order, where the 

two different symmetry breakings are switched. 

Let us discuss the nature of the spin-imbalance phase 
in more detail. As we disucssed above, the imbalanced 
order is uniform along the leg direction, while the field 
theory based on the effective model @ suggests the emer- 
gence of a staggered imbalance order ((/Xj) = —{fjLj+i)). 
This mismatch of the effective theory may be attributed 
to the fact that the imbalanced order is located at very 
vicinity of the ferro-chirality transition point J^2i where 
the velocity v almost vanishes and thus the system be- 
comes fragile. Furthermore, we find that the rapid in- 
crease of jj, in J' > J^^ causes a rapid raise of the en- 
ergy of the unit triangle (DMRG data are not presented 
here), implying that the effect of J' nonperturbatively 
reduces the energy of the intra-triangle bonds. Thus it 
is suggested that the role of the intra-triancle coupling 
becomes essential and thus the T — ^ sector certainly 
hybridizes into the plateau state in the spin-imbalance 
phase. 

The effective model is based on the massive weight 
of the T = ^ sector, while the mixing of the T = | sec- 
tor is possibly essential for the spin imbalance phase. We 
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FIG. 4. (color online) Spin profile (Sf j) of the spin imbalance 
phase. The tube length is L = 96 and the inter triangle 
coupling is J' — 0.485. Solid circles denote the expectation 
value of the most polarized spin in each rung triangle and 
the open circles correspond to those of remainsing two spins 
on the triangle. The horizontal broken lines is the averaged 
magnetization of each rung in the plateau state. 
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FIG. 5. (color online) Correlation function {{SfjSi^,) — | 
in the spin imbalance phase. Solid circles is the correlator for 
the chain consisting of the most polarized spins on the unit 
triangules, and the open circles correspond to that for the 
remaining two chains. Inset represents J' dependence of the 
inverse correlation length along the chains for the most 
and less polarized spins in the triangle. 
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should thus investigate the expectation value of — 
(T^/3 — 5/4), which is the projection operator into the 



1/2, _ 



1 and it 



T = i sector. In the J' = limit, {P. 

gradually decreases up to J'^^. Figure H] shows (P^/^) 
around the transition points, which is obtained by the 
infinite DMRG. In the figure, we can see that the be- 
havior of (P^/^) drastically changes at Jd ~ 0.478 and 
Jc2 = 0.496. In Jci < J' < Jc2, (P^/^) rapidly decreases 
with increasing J'. This supports that the driving mech- 
anism of the spin-imbalance phase relies on the mixing 
of the T — ^ sector. Although a Berezinskii-Kosterlitz- 
Thouless (BKT) type transition accompanying the Z3 
symmetry breaking^^ is naively expected at Jd, the na- 



FIG. 6. (color online) Expectation value of the projection 
operator (P^^^) on the unit triangle. 



ture of the phase transition might be essentially modi- 
fied by the T — ^ sector. However, we may claim within 

the present analysis that(P^/^) is continuously changed 
around Jd , suggesting a continuous quantum phase tran- 
sition. Further analysis is necessary to completely de- 
termine the nature of this transition, including the uni- 
versality class. On the other hand, there exists a clear 
jump of (P^/^) at Jc2 = 0.496. The two branches near 
Jc2 represent two self-consistent solutions corresponding 
to the chirality ordered and spin-imbalance states in the 
DMRG iterations; the solution of the previous parameter 
is used as an initial state for the next parameter, so that 
the metastable states can be reproduced. By comparing 
energies of the two branches, the first-order transition 
point can be determined as — 0.496. This result is 
consistent with the behaviors of the order parameters in 
Fig. [31 



IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have explored the quantum phase 
transitions of the ^ plateau state of the spin tube. In con- 
trast to the usual plateaus of one-dimensional spin sys- 
tems (chains and ladders) , the chirality degree of freedom 
generated from the tube structure plays crucial roles. 
The results are summarized in Fig. [T] (b), where the 
chirality liquid phase with gapless non-magnetic excita- 
tions, the spin-imbalance phase and the ferro-chirality 
phase emerge. The qualitative features of these phases 
may be explained by the effective chirality model (U) and 
the S'3-symmetry breakings. However, the precise analy- 
sis of the projection operator P^/^ has revealed that the 
uniform spin imbalance order is driven by mixing of the 
T = I sector, which is beyond the scope of the effective 
model (|4]). The transition between the chirality liquid 
and the spin-imbalance phase is of continuous type, and 
the fluctuation of less polarized spins in the imbalance 
phase becomes divergent near the transition. On the 



6 



other hand, the transition between the spin-imbalance 
and ferro-chirahty ordered phases is shown to be of first 
order type. 

Here it should be commented that another spin- 
imbalance phase with gapless magnetic excitations is ex- 
pected in a high magnetic fieldi^. Its connection to 
the present spin-imbalance phase may be an interest- 
ing problem for through understanding of mechanisms 
of the spin imbalance. As we mentioned in the introduc- 
tion, a chirality-ordered spin liquid appears in the weak 
rung-coupling region J ^ J' in magnetic fields^^. This 
spin liquid is expected to change into the 1/3 plateau 
state with the chirality order—' via a BKT transitionii at 
the order of J / J' — 0.1. Combining our present results 
with this, we can conclude that, as J' increases from the 
strong rung limit, the chirality liquid, spin-imbalance or- 
der, ferro-chirality order, and ferro-chirality-ordered spin 
liquid can be observed at m = ^ in order. 

An important aspect of the spin tube is that the phase 
transitions occur without destroying the plateau. The 
energy scale of the chirality is significantly lower than 
the width of the large plateau. Therefore, for exam- 
ple, a specific heat measurement will solely observe a 
linear temperature dependence originating from the chi- 
rality modes in the wide spin-gapped plateau region of 
J' < J^]^ , in contrast to the twisted tube^^. From experi- 



mental viewpoint, moreover, another plausible feature of 
the spin tube is that the gapped chirality order is ex- 
panded in the wide range of J', which is contrasted to 
the narrow chirality-ordered phases with gapped mag- 
netic excitations in the classical XY model on triangu- 
lar lattice- and spin-S* J1-J2 chain o^^'^^ . If a coupling 
between chirality and electric polarization is introduced, 
the chirality order can induce a ferro-electric polarization 
in spite of the absence of any magnetic ordering. Also, 
the similar chirality degree of freedom is discussed in the 
coupled trimer model, which may reduce to the spin tube 
in an anisotropic limit We thus believe that the spin 
tube provides a fascinating play ground of the chirality 
degrees in the realistic experimental situation. 
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